## make a grid
## use the ice grid
library(raster)
## Loading required package: sp
library(aceecostats)
#grid <- raster(extent(-3950000, 3950000, -3950000, 4350000), nrow = 332, ncol = 316,
# crs = "+proj=stere +lat_0=-90 +lat_ts=-70 +lon_0=0 +k=1 +x_0=0 +y_0=0 +a=6378273 +b=6356889.449 +units=m +no_defs")
library(dplyr)
##
## Attaching package: 'dplyr'
## The following objects are masked from 'package:raster':
##
## intersect, select, union
## The following objects are masked from 'package:stats':
##
## filter, lag
## The following objects are masked from 'package:base':
##
## intersect, setdiff, setequal, union
grid <- extent(aes_zone) %>% tabularaster::bufext(25000) %>% raster(res = 25000, crs = projection(aes_zone))
library(feather)
library(sf)
## Linking to GEOS 3.5.0, GDAL 2.1.2, proj.4 4.9.2
data("wrld_simpl", package= "maptools")
m <- st_geometry(st_transform(st_as_sf(wrld_simpl), st_crs(projection(grid))))
outf <- "/mnt/acebulk"
vars <- c("sst", "chl")[1]
for (ivar in seq_along(vars)) {
obj <- brick(file.path(outf, sprintf("%s.grd", vars[ivar])))
#obj <- projectRaster(obj, grid, method = "bilinear")
obj <- setZ(obj, getZ(brick(file.path(outf, sprintf("%s.grd", vars[ivar])))))
segs <- cumsum(c(0, abs(diff(unclass(factor(aes_season(getZ(obj))))))))
cell_tab <- vector("list", length(unique(segs)))
dates <- as.POSIXct(getZ(obj))
for (i in seq_along(cell_tab)) {
asub <- which(segs == unique(segs)[i])
a_obj <- setZ(readAll(subset(obj, asub)), dates[asub])
print(min(dates[asub]))
gridmax <- projectRaster(max(a_obj, na.rm = TRUE), grid)
gridmean <- projectRaster(mean(a_obj, na.rm = TRUE), grid)
gridcount <- projectRaster(calc(a_obj > 0, sum, na.rm = TRUE), grid, method = "ngb")
plot(as(extent(gridmax), "SpatialPolygons"), col = "black", asp = 1); plot(gridcount, col = viridis::viridis(100), add = TRUE)
print(i)
rm(tab, a_obj)
gc()
}
}
## [1] "1981-09-01 UTC"
## [1] 1
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1981-12-01 UTC"
## [1] 2
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1982-03-01 UTC"
## [1] 3
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1982-06-01 UTC"
## [1] 4
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1982-09-01 UTC"
## [1] 5
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1982-12-01 UTC"
## [1] 6
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1983-03-01 UTC"
## [1] 7
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1983-06-01 UTC"
## [1] 8
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1983-09-01 UTC"
## [1] 9
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1983-12-01 UTC"
## [1] 10
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1984-03-01 UTC"
## [1] 11
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1984-06-01 UTC"
## [1] 12
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1984-09-01 UTC"
## [1] 13
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1984-12-01 UTC"
## [1] 14
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1985-03-01 UTC"
## [1] 15
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1985-06-01 UTC"
## [1] 16
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1985-09-01 UTC"
## [1] 17
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1985-12-01 UTC"
## [1] 18
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1986-03-01 UTC"
## [1] 19
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1986-06-01 UTC"
## [1] 20
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1986-09-01 UTC"
## [1] 21
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1986-12-01 UTC"
## [1] 22
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1987-03-01 UTC"
## [1] 23
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1987-06-01 UTC"
## [1] 24
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1987-09-01 UTC"
## [1] 25
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1987-12-01 UTC"
## [1] 26
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1988-03-01 UTC"
## [1] 27
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1988-06-01 UTC"
## [1] 28
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1988-09-01 UTC"
## [1] 29
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1988-12-01 UTC"
## [1] 30
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1989-03-01 UTC"
## [1] 31
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1989-06-01 UTC"
## [1] 32
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1989-09-01 UTC"
## [1] 33
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1989-12-01 UTC"
## [1] 34
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1990-03-01 UTC"
## [1] 35
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1990-06-01 UTC"
## [1] 36
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1990-09-01 UTC"
## [1] 37
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1990-12-01 UTC"
## [1] 38
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1991-03-01 UTC"
## [1] 39
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1991-06-01 UTC"
## [1] 40
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1991-09-01 UTC"
## [1] 41
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1991-12-01 UTC"
## [1] 42
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1992-03-01 UTC"
## [1] 43
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1992-06-01 UTC"
## [1] 44
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1992-09-01 UTC"
## [1] 45
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1992-12-01 UTC"
## [1] 46
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1993-03-01 UTC"
## [1] 47
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1993-06-01 UTC"
## [1] 48
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1993-09-01 UTC"
## [1] 49
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1993-12-01 UTC"
## [1] 50
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1994-03-01 UTC"
## [1] 51
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1994-06-01 UTC"
## [1] 52
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1994-09-01 UTC"
## [1] 53
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1994-12-01 UTC"
## [1] 54
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1995-03-01 UTC"
## [1] 55
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1995-06-01 UTC"
## [1] 56
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1995-09-01 UTC"
## [1] 57
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1995-12-01 UTC"
## [1] 58
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1996-03-01 UTC"
## [1] 59
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1996-06-01 UTC"
## [1] 60
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1996-09-01 UTC"
## [1] 61
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1996-12-01 UTC"
## [1] 62
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1997-03-01 UTC"
## [1] 63
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1997-06-01 UTC"
## [1] 64
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1997-09-01 UTC"
## [1] 65
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1997-12-01 UTC"
## [1] 66
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1998-03-01 UTC"
## [1] 67
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1998-06-01 UTC"
## [1] 68
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1998-09-01 UTC"
## [1] 69
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1998-12-01 UTC"
## [1] 70
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1999-03-01 UTC"
## [1] 71
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1999-06-01 UTC"
## [1] 72
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1999-09-01 UTC"
## [1] 73
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "1999-12-01 UTC"
## [1] 74
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2000-03-01 UTC"
## [1] 75
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2000-06-01 UTC"
## [1] 76
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2000-09-01 UTC"
## [1] 77
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2000-12-01 UTC"
## [1] 78
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2001-03-01 UTC"
## [1] 79
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2001-06-01 UTC"
## [1] 80
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2001-09-01 UTC"
## [1] 81
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2001-12-01 UTC"
## [1] 82
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2002-03-01 UTC"
## [1] 83
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2002-06-01 UTC"
## [1] 84
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2002-09-01 UTC"
## [1] 85
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2002-12-01 UTC"
## [1] 86
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2003-03-01 UTC"
## [1] 87
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2003-06-01 UTC"
## [1] 88
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2003-09-01 UTC"
## [1] 89
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2003-12-01 UTC"
## [1] 90
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2004-03-01 UTC"
## [1] 91
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2004-06-01 UTC"
## [1] 92
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2004-09-01 UTC"
## [1] 93
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2004-12-01 UTC"
## [1] 94
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2005-03-01 UTC"
## [1] 95
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2005-06-01 UTC"
## [1] 96
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2005-09-01 UTC"
## [1] 97
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2005-12-01 UTC"
## [1] 98
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2006-03-01 UTC"
## [1] 99
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2006-06-01 UTC"
## [1] 100
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2006-09-01 UTC"
## [1] 101
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2006-12-01 UTC"
## [1] 102
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2007-03-01 UTC"
## [1] 103
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2007-06-01 UTC"
## [1] 104
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2007-09-01 UTC"
## [1] 105
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2007-12-01 UTC"
## [1] 106
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2008-03-01 UTC"
## [1] 107
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2008-06-01 UTC"
## [1] 108
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2008-09-01 UTC"
## [1] 109
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2008-12-01 UTC"
## [1] 110
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2009-03-01 UTC"
## [1] 111
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2009-06-01 UTC"
## [1] 112
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2009-09-01 UTC"
## [1] 113
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2009-12-01 UTC"
## [1] 114
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2010-03-01 UTC"
## [1] 115
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2010-06-01 UTC"
## [1] 116
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2010-09-01 UTC"
## [1] 117
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2010-12-01 UTC"
## [1] 118
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2011-03-01 UTC"
## [1] 119
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2011-06-01 UTC"
## [1] 120
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2011-09-01 UTC"
## [1] 121
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2011-12-01 UTC"
## [1] 122
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2012-03-01 UTC"
## [1] 123
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2012-06-01 UTC"
## [1] 124
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2012-09-01 UTC"
## [1] 125
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2012-12-01 UTC"
## [1] 126
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2013-03-01 UTC"
## [1] 127
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2013-06-01 UTC"
## [1] 128
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2013-09-01 UTC"
## [1] 129
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2013-12-01 UTC"
## [1] 130
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2014-03-01 UTC"
## [1] 131
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2014-06-01 UTC"
## [1] 132
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2014-09-01 UTC"
## [1] 133
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2014-12-01 UTC"
## [1] 134
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2015-03-01 UTC"
## [1] 135
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2015-06-01 UTC"
## [1] 136
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2015-09-01 UTC"
## [1] 137
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2015-12-01 UTC"
## [1] 138
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2016-03-01 UTC"
## [1] 139
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2016-06-01 UTC"
## [1] 140
## Warning in rm(tab, a_obj): object 'tab' not found

## [1] "2016-09-01 UTC"
## [1] 141
## Warning in rm(tab, a_obj): object 'tab' not found
